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In this paper we examine the cosmological consequences of fourth order Galileon gravity. We 
carry out detailed investigations of the underlying dynamics and demonstrate the stability of one de 
Sitter phase. The stable de Sitter phase contains a Galileon field tt which is an increasing function 
of time (7i > 0). Using the required suppression of the fifth force, supernovae, BAO, and CMB 
data, we constrain parameters of the model. We find that the tt matter coupling parameter f3 is 
constrained to small numerical values such that f3 < 0.02. We also show that the parameters of the 
third and fourth order in the action (C3, C4) are not independent and with reasonable assumptions, 
we obtain constraints on them. We investigate the growth history of the model and find that the 
sub-horizon approximation is not allowed for this model. We demonstrate strong scale dependence 
of linear perturbations in the fourth order Galileon gravity. 

PACS numbers: 



I. INTRODUCTION 

According to the dark energy paradigm, the phe- 
nomenon of late time cosmic acceleration |l| can be un- 
derstood by assuming that an exotic relativistic fluid with 
large negative pressure fills the whole Universe [2]. The 
simplest example of such a homogeneous fluid is pro- 
vided by cosmological constant A which is automatically 
present in the Einstein equations by virtue of Bianchi 
identities. The dark energy model based upon the cos- 
mological constant a la AC DM is consistent with all the 
cosmological observations at present. However, the field 
theoretic understanding of A is far from being satisfac- 
tory and its small numerical value leads to the well-known 
coincidence and fine tuning problems. A variety of scalar 
fields as candidates of dark energy were then investigated 
in the hope of addressing the said problems. The scalar 
field models with generic features might alleviate the fine 
tuning and coincidence problems leaving the dark energy 
metamorphosis as a challenge for future observations in 
cosmology. 

There is an alternative view which advocates the neces- 
sity for the paradigm shift according to which late time 
cosmic acceleration is an artifact of large scale modifi- 
cation of gravity rather than the consequence of dark 
energy. On theoretical grounds it is plausible that grav- 
ity suffers modifications at large scales where it is never 
tested directly. We know that gravity gets corrected 
quantum mechanically at small scales which is beyond 
the observational reach at present. As for the large scale 
modification, it should give rise to late time acceleration, 
be distinguishable from A and at the same time be consis- 
tent with local gravity constraints. The last requirement 
is nontrivial as Einstein gravity agrees with solar physics 
and the equivalence principal with a high accuracy. 

There are broadly two ways used to evade the local 
gravity constraints namely the chameleon mechanism Q 
and Vainstein screening Q. The first method is widely 
used in f(R) theories of gravity [E| with a disappearing 



cosmological constant. In this scenario, the mass of the 
scalar degree of freedom dubbed scalaron present in the 
theory becomes large thereby hiding the scalaron locally. 

In generic models of f(R) gravity [6], the chameleon 
mechanism allows us to satisfy the local gravity con- 
straints but at the same time makes these models vul- 
nerable to the problem of a curvature singularity whose 
resolution requires the fine tuning worse than the one 
encountered in ACDM model. The problem can be alle- 
viated by invoking an R 2 correction but the scenario be- 
comes problematic if extended to the early universe. The 
later puts an stringent constraint on a viable large scale 
modification of gravity within the framework of f(R) the- 
ories. 

The scalar degrees of freedom also naturally arise in 
the four dimensional effective theories. They couple to 
matter source and might give rise to fifth force effects. 
In this case, the local gravity constraints can be evaded 
using the Vainstein screening mechanism which is im- 
plemented by using the nonlinear self interaction of the 
scalar field. Nonlinearity becomes important in the vicin- 
ity of dense objects allowing the scalar degrees to decou- 
ple from the matter source. In DGP [7], the scalar degree 
of freedom appears in the form of brane bending mode 
with the required nonlinear derivative interactions of the 
simplest type which is invariant under the shift symme- 
try in the flat space time. The equation of motion for 
the scalar field dubbed Galileon is necessarily of second 
order. The general structure of higher order Lagrangian 
of the Galileon field was obtained in [8] . Similar to Love- 
lock gravity, the Galileon gravity provides a consistent 
modification of GR leaving the local physics intact. The 
DGP model which includes the lowest order Galileon La- 
grangian in its decoupling limit suffers from the problem 
of instabilities [9]. The Galileon modified gravity in its 
general setting can give rise to late time acceleration and 
is free from negative energy instabilities. The model has 
the well-posed Cauchy problem and is safe from para- 
doxes related to micro-causality [lOj . 
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In this paper we investigate the cosmological dynam- 
ics based upon Galileon gravity, set up the autonomous 
system and discuss the existence and stability of fixed 
points. We especially focus on the self-accelerating so- 
lution and explore the observational constraints on the 
model parameters using supernovae, BAO and CMB 
data. We also study metric perturbations and investi- 
gate the growth history of the model. 

II. BACKGROUND 

In Galileon theories, the large scale modification of 
gravity arises due to the nonlinear derivative self interac- 
tion of a scalar field tt dubbed the Galileon field, which 
couples with matter and metric. In what follows, we shall 
consider the Galileon action of the form, 




where {q} are constants and the are given by 



With A = C 3 ~ 12C 2 C 4 ±C 3 y/ 4 -8C 2 C 4 

=•= c 4 

We have, therefore, two de Sitter solutions for this 
model, namely, the positive branch (A+) and the 
negative branch (A-). 

In Ref. [To| , various conditions of stability of the 
theory were derived. It was shown that positive values 
for the parameters (c2,C4,/3) and cs > ^80204 can give 
rise to viable evolution. 

It is straightforward to show that A- < 0. As for the 
A_, it is a decreasing function of C3 and therefore the 
highest value of A- is achieved when C3 = \j8c2C4. We 
found A™ ax = -4c 2 < 0. 

If we consider the conditions of the stability of the 
theory, then the negative branch is ruled out as it does 
not have a de Sitter phase thereby leaving us with one 
self-accelerating solution in the positive branch. Bearing 
this in mind, in the discussion to follow, we shall only 
consider the positive branch of the theory and thus 
redefine A = A+. 
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-i(V7r) 2 D7r 
-i(V7r) 2 [(Dtt) 2 - Tr^Tr"" + nVG^] 
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The fourth order Galileon theory leads to the following 
evolution equations in the FRW background [Iq| , 



III. AUTONOMOUS SYSTEM 

For the sake of convenience, we use the system of units 
such that 7r is dimensionless and so is also true for C2 . 
We then define two new dimensionless parameters £3 = 
csH^^ C4 = c^H^g where H^s is the Hubble function 
during the de Sitter era. It can easily be noticed that the 
parameters (C2, £3, £4) are not independent. In fact if we 
consider the de Sitter solution derived in Ref. [lOj, we 
find that, 



3H 2 



p r + — tt 2 - 3c 3 Htt 6 + —c A H' z ^ (6) 
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2H + SH 2 



1 c 2 . 2 . 2 » 

-3/V - —7T -C 3 7T 7T 

-^c 4 tt 3 (3H 2 7T + 2HTT + 8#7f) (7) 



Ppm = -C 2 (SHJT + TT) + 3C 3 7T (3H 2 7T + HjT + 2#7T 

-18C4#7T 2 ^ 2 7T + 2#7T + 3#7f) , (8) 

where H = a/a is the Hubble function, p m and p r are 
the density of matter and radiation respectively 

It was found [To[ that the model has a self-accelerating 
solution iff 



c| - 8c 2 c 4 > 0, 
A+ > or A_ > 



(9) 
(10) 



C 3 + 8c 2C 4 
H-dS^dS = — 
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It is then straightforward to find the relation between 

(c 2 ,c 3 ,c 4 ), 



48 
~A 



^3 + \Zcj ~ 8c 2 c 4 
12c 4 



(13) 



It is therefore obvious that the parameters are not in- 
dependent if we normalized them by H^s- On the other 
hand we can define a normalization of these parameters 
at any redshift, in particular we can do it today by 



rr-2 ( Hds\ , 
C3 = C 3 H =[ — 1 C 3 



c 4 c 4 # 
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The factor H^s/Ho depends on the evolution of the 
model but for an evolution close to ACDM, we have 



c 3 
c 4 



Q A c 3 



(16) 
(17) 



For physical interpretation we will work with (£3,(34) 
given by expressions (|14I15|) . In this case, it is reasonable 
to assume that the parameters (£3,04) are of order one. 
In fact it was showed in [loj that the Galileon force is 

suppressed for scales smaller than r* where r * = \/ ftr 



c 4 

2/3 Hi 

and r s is the Schwarzschild radius of the source. If we 
consider C4 and j3 of the same order, we have 



r* (Earth) = 1.7 pc 



r*(Sun) = 120 pc 



(18) 



(19) 



TV (Milky Way) = 1.2 Mpc (20) 

which are the same scales as in the DGP model. 

The evolution equations (|6I7|8|) can easily be cast in 
the autonomous form. 

Let x = tt/H and y = ttH/Hq with Hq as the Hubble 
constant today. The evolution equations acquire the form 
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where a prime represents a derivative with respect to 
N = In a and 



7T 3y c 2 x (4 + x 2 y(c 3 - 18c 4 y)) - (2 - 3c 4 x 2 y 2 ) (3c 3 xy - 2pVt m ) + 4xy(c 3 - 12c 4 y)Q r 
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2x 54:C4X 2 y 3 (c 3 — 5cav) + c-2 (—2 + 3c±x 2 y 2 ) — 3y (3§c±y + c 3 (—4 + c 3 x 2 y)) 

Ax,y 

108c4X 2 ?/ 3 (c3 — 5c42/) + C2 (—4 + 6c4X 2 y 2 ) — 6y (3§c/±y + C3 (—4 + c 3 x 2 y)) 

1 - tt r - \x 2 (c 2 - 6c 3 y + 45c 4 2/ 2 ) 
6 v 7 
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3x 2 7/(-4c 3 



6c 32 / (6 - 
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45c 4 £ ?/ 



39c 4 ?/) + 4Q r ] 
xf3Vt 
4xPQ m + 12^ r 



45c 4 x 2 ?/ 2 



4ft r ) 
(27) 



Given the definition of the two variables (x,y), it 
is enough to consider the phase space subject to the 
condition, sign(x) = sign(y). 

This autonomous system has a saddle point, which cor- 
responds to the radiation era: 



P ril :(x,y,n r ) = (0,0,1) 



(28) 



The fixed point which corresponds to the matter- 
dominated epoch is 



P m : (x,y,Q r ) = (0,oo,0), (29) 
whereas the two de Sitter points are given by 



P dS : (x,y,Q r ) 



/48 c 3 + yci - 8c 2 c 4 
1 V A ' 12c 4 



0) (30) 



We note that yds is always positive because of the con- 
ditions of stability of the theory (£4 > 0,C3 > V8C2C4). 
Hence, we shall consider only those solutions which sat- 
isfy, > (sign(x dS ) = sign(y dS )). In this case, the 
system has only one de Sitter point given by, 



8C264 



,0) (31) 



which is an attractor iff A > 16/3 2 /3. In the case 
of physical interest with a small numerical value of 
f3 , we recover the condition (A > 0) derived in [loj . 
In this last case, the de Sitter point is always an attractor. 

In addition to these three eras, various other critical 
points can be found. The relevant points are 
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P[ : (x, 2/, ft r ) = (±y — >0,0), = 0, w e ff = 1, saddle point 

02 



Pi '• (x,y,Q r ) = ( ?0,0), ^eff = j = 1 , attractor if f3 2 < — and a saddle point otherwise. 

C2 3C2 3C2 2 



P 3 : (x,2/,O r ) = (-, W — — — ,0) , ^eff 



-l, n ri 



-4-9 



c 2 c 3 /9c 2 + 6/3 2 



, which is a saddle point 



p 2 /3 2 y 

or an attractor depending on the set of parameters (c2, C3, £4, /?). 
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In Figs. [U and [2j we show the evolution of the au- 
tonomous system. We have a standard evolution in the 
first case where we chose positive initial conditions for x 
and y. The system evolves along the axis x = during 
the matter phase before it is attracted by the de Sitter 
point along the line y = yds- We do not have a track- 
ing solution when f3 7^ 0, but after the matter phase the 
coupling to matter j3 is weak therefore we recover the 
tracking solution y = yds derived in 
In the second case we chose negative initial conditions, in 
this case there is no de Sitter point in the subspace con- 
sidered, after a matter phase the model is attracted by 
P2 (the attractor for these values of the parameters). We 
can see in the bottom plot of Fig. [2] that in the future, 
Q m = 1 — 2^- ~ 2/3 which corresponds to the attractor 
P 2 - 

It is clear that we can not achieve a de Sitter phase if 
we consider initial conditions in the subspace (x, y) < 0. 
Therefore, we have to consider positive initial conditions, 
which means tt > 0. 



IV. OBSERVATIONAL CONSTRAINTS 

As demonstrated in the previous section, the param- 
eters of the model are not free; we can consider 04 as a 
function of (/3, 02,03). Thus in this section, we shall im- 
pose the constraints on the parameters (/?, 02,03) only. 
We constrain the parameters of the model by using su- 
pernovae, BAO and CMB datas. We used the compiled 
constitution set of 397 type la supernovae [12| for which 
the x 2 is defined as 
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FIG. 1: Top Panel: The projected phase space in the plane 
(x,y) in Poincare coordinates for f3 = 0.1, C2 = 1, £3 = 15, 
C4 = 4. The circles represent critical points. The initial con- 
ditions are chosen in the subspace (x, y) £ K+. Bottom Panel: 
The evolution of Q, as a function of log(l + z). 



2 _ (/fthj ~ Mobs,i) 2 

XsNla — ~2 



(32) 



with 



where /io = 25 + log H 1 is marginalized [HI, [HI and 
d,L is the luminosity distance. 



Vth,i = Slog d L (zi) +/i 



We used the BAO distance ratio D v (z = 0.35)/ D v (z 
(33) 0.2) = 1.736 ± 0.065 0, where 
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FIG. 2: Top Panel: The projected phase space in the plane 
(x,y) in Poincare coordinates for f3 = 0.7, c<i = 1, C3 = 15, 
C4 = 4. The circles represent critical points. The initial condi- 
tions are chosen in the subspace (x,y) E M_. Bottom Panel: 
The evolution of O as a function of log(l + z). 



D v {z) 



H(z) 



dz' 



1/3 



(34) 



Finally we used the CMB shift parameter R = 1.725 ± 
0.018 [14J , where 



R = H { 



JO 



dz' 
H{z) 



(35) 



We observe (see Fig. (j3])) that when C2 = /3 (motivated 
by a Galileon force of the order of the gravitational force 
at large scales [To| ) , /3 is constrained by the data to small 
values j3 < 0.02, however the parameter cs is constrained 
at 2<j at small values but unconstrained at 3a. For the 
model where we fixed C2 = 1 we found that parameters 
are constrained to take small values even at 3a. 
In the bottom plot of the same figure, we chose a small 
value for the coupling f3 = 0.01, we observe that for a 
fixed C2 large values of cs are preferred. 





- 




- 


: V 1 






FIG. 3: Top Panel: Contour plots at la, 2a and 3a for C2 = /3. 
Middle Panel: Contour plots at la, 2a and 3a for C2 = 1. 
Bottom PanehContour plots at la, 2a and 3a for f3 = 0.01. 
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V. PERTURBATIONS 

As mentioned earlier, the model is constrained differ- 
ently by the data, when the assumption on c 2 is differ- 
ent. But we can always find a range of parameters which 
can fit the data, this is obvious because the model has a 
matter phase and a self-accelerating solution as for the 
ACDM. Therefore the large number of free parameters 
compared to the ACDM model cannot be constrained. 
It is widely discussed in the literature that we can have 
a strong signature of modified gravity models by looking 
to the evolution of perturbations [18[. We performed 
this analyze in the linear regime. We will consider, 
k < 0.1/iMpc -1 , in this case we are at scales larger than 
r*, where the linear approximation is allowed. 
We will consider the perturbations in the comoving 
gauge. Therefore, we can write the equations of per- 
turbations for (5 m ,<$7r) as the following: 



6m + AiS m + A 2 p m S m + AsSn + A^Stt + A 5 Stt 

= (36) 

B x 'S m + B 2 S m + BsPmSm + B 4 Stt + B 5 Stt + B 6 6tt 

= (37) 

The coefficients {Ai,B{} are given in the Appendix. 
In the subhorizon approximation, we can approximate 

Eq. (jail) b y 



A 2 Bi — B3 
H = B e - Bl A 5 pmS " 



(38) 



and this gives the equation for S m under this approxi- 
mation 



A^ 



(39) 



with 



A 2 B\ — £?3 

Bq — B\Ar, 



(40) 



During the matter phase x = n/H <C 1 and y — 
ttH/Hq ^> 1, therefore we have 



B 6 



2H, 
-1/2, 



-If » 1 

Then Eq. (|39|) can be approximated by 



2H5„ 



0, 



(41) 
(42) 
(43) 



(44) 
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FIG. 4: Geff versus the redshift log(l+z) for f3 = 0.01, c 2 = 1, 
C3 — 15 and C4 — 4. 



We recover the equation of perturbations of standard 
general relativity during the matter phase. The evolu- 
tion of S m for the Galileon model is equivalent to the 
evolution of matter perturbation in the standard model 
of cosmology. However after the matter phase, the 
approximation x = tt/H <C 1 and y = ttH/Hq ^> 1 is 
no longer valid and a deviation from the ACDM model 
appears. 

In Fig. [Hwe can see that G e ff — 1 until z ^ 3.5. But 
at small redshift G e fr diverges and the asymptotic value 
is negative. It does not necessarily imply the instability 
of the theory; this may be due to the subhorizon approx- 
imation that we did more than a problem of the theory. 
It was shown in Ref. [17] that replacing a second-order 
differential equation by an algebraic equation (f37j) can be 
misleading. 

Therefore we have to solve the complete system of 
coupled differential equations (j36j37j). 

Figure [5] shows 7 as a function of the redshift for the 
following parameters (/?, c 2l £3, £4) = (0.01, 1, 15, 4) where 
7 is defined by 



7 



din S m 
din a 



(45) 



For all the modes considered in the linear regime 
(k < 0.1/iMpc _1 ), 7 has an oscillating mode during 
the matter phase; this oscillation becomes negligible for 
z < 1.5. We remark that 7 is crucially different from its 
counterpart in models of dark energy within the frame- 
work of general relativity or modified gravity models 
as f(R) or chameleon gravity. We have a strong dis- 
persion of 70 with 70 (fc = 0.1 h Mpc -1 ) = —4.4 and 
7o(/c = 0.001 h Mpc" 1 ) = 0.19. 
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FIG. 5: 7 versus log(l + z) for f3 = 0.01, c 2 = 1, c 3 = 15 and 
C4 — 4. 
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FIG. 6: 70 = j(z = 0) for different values of (02,03) and 
P = 0.01 for the scale k = 0.001 h Mpc -1 . 



(^2? C3), 70 is shifted by —4 if we consider a smaller scale 
k = 0.1 h Mpc" 1 . 



VI. CONCLUSION 

In this paper, we have shown that the Galileon modi- 
fied gravity has only one stable de Sitter branch dubbed 
positive branch. We find that tt should be positive for 
the de Sitter phase to be reached. We have analyzed the 
stable branch by considering the constraints from Super- 
novae, BAO and CMB and found that the coupling pa- 
rameter f3 is constrained to small values, f3 < 0.02. The 
parameters (C2, £3, £4) are not independent. Therefore we 
restricted our analysis to two parameters namely, (C2, £3) 
. We found that the contour plots are strongly dependent 
on the assumption. For C2 = 1, motivated by a standard 
normalization of the kinetic term in the theory , we found 
that the model is strongly constrained even at 3a. We 
also repeated our calculations for C2 = /3 which is moti- 
vated by the requirement that the Galileon force should 
be of the order of the gravitational force at large scales. 
For this last case, the model is constrained at a small 
value of the parameters at 2a but C3 and therefore, C4 is 
unconstrained at 3a. Finally, we performed the analysis 
of linear perturbations in the model under consideration. 
We found that the subhorizon approximation can not be 
considered for the Galileon model, we have to consider 
the perturbations of the field too. The growth of linear 
perturbations 7 turns out to be strongly dependent on 
the scales considered. At very large scales 70 — 0.25 and 
this value is negative for smaller scales which can rule out 
the model. It should be noted that the addition of the 
fifth term in the Galileon theory can produce different 
results as it contains an additional free parameter (05). 
The inclusion of field potential can also be considered to 
avoid the unwanted behavior in the linear regime of per- 
turbations. In this case, the model can mimic features 
similar to the scenarios with a chameleon mechanism. In 
our opinion, it might be interesting to investigate this 
possibility. 



In the case considered in Fig. [6] we showed the varia- 
tion of 70 = j(z = 0) for < c 2 < 2 and < c 3 < 10. We 
find that at small values of C2, 70 is going larger (70 > 0.2) 
for the scale k = 0.001 h Mpc -1 . In the same range of 
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Appendix A: Perturbations 
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Ai = 2 [ 8H + ^ 3 (- 2c 3 + c 4 (-8il7r + 42c 4 #tt 5 - 24tt - 7T 4 (c 3 - 36c 4 7t)))] (Al) 

(2 — 3c 4 7T 4 ) 

A 2 = - 2 + C4 ^ 2 (A2) 

(2-3C47T 4 ) 2 V 

A 3 = -/3 + 37T _ 2 + 3c4 . 4 (A3) 
A A = -3H/3 + 2 [-8i?/3 + 80c 4j ff/37r 4 - 150c^iJ/37r 8 - 4tt (2c 2 + 3c 3 7f) + 2c 4 tt 5 (4c 2 + 9c 3 7f) 

(2 - 3c 4 7T 4 ) 

+9c 4 iJ7r 6 (c 3 - 36c 4 7f) + 18Htt 2 (c 3 + 12c 4 tt) - 4tt 3 (l8c 4 if 2 + c 3/ 3 - 18c 4 H - 6c 4 /?tt) 

-2c 4 7r 7 (l26c 4 if 2 - 5c 3/ 3 + 54c 4 ij + 18c 4/ 3tt) (A4) 



M= 1 



(2 - 3c 4 tt 4 ) 2 



A: 2 

-2 (4/3 - 16c 4 #7T 3 - 12c 4 /3tt 4 - 24cIHtt 7 + 9c^tt 8 + c 4 tt 6 (c 3 - 36c 4 tt) + 2tt 2 (c 3 + 12c 4 7f)) 



-/3 (24H 2 + 48c 4 # 2 tt 4 + 378cjH 2 7r 8 + 3# (4 - 36c 4 tt 4 + 45c^tt 8 ) + 8tt 2 (c 2 + 3c 3 7f ) - 4c 4 tt 6 (2c 2 + 9c 3 tt) 
-12c 4 #tt 7 (c 3 - 54c 4 7f) - 24#tt 3 (c 3 + 18c 4 7r))] - 3/3# (A5) 

Si = -C 3 7T 2 + 12c 4 #7T 3 (A6) 

/ / o • \ 9 , N \ 4c 4 tt 2 (-2H - c 3 tt 3 + 15c 4 #tt 4 ) (#tt + 3tt) , 
£ 2 = tt (c 2 + 6c 4 (9tf 2 + 2tf J ^ 2 - 2c 3 7f - 6Htt (c 3 - 6c 4 7f)) = * _ 2 - ,T 4 J - (A7) 

B 4 = -c 2 + QcsHrr + (c 3 /3 - 54c 4 # 2 ) tt 2 - 12c 4 #/3tt 3 (A9) 

4c 4 tt 2 1-2HP + c 2 tt - 9c 3 #tt 2 + (90c 4 # 2 - c 3 /3) tt 3 + 15c 4 #/3tt 4 ) (Htt + 3tt) 
£5 = -2 + 3C47T* 3S(c 2 -2c 3 if) 

-3#tt 2 (54c 4 # 2 - 5c 3 /3 + 36c 4 # + 12c 4 f3jr^j + 2tt (9c 3 # 2 - c 2 /3 + 3c 3 # + (c 3 /3 - b4c 4 H 2 ) tt) 

-12c 4 /3 (l2# 2 + ij) tt 3 (A10) 

2 f 26c 4 # 2 + c 3 /3 + 12c 4 #) tt 2 - 24c 4 Hp7r 3 - 3c 4 f 10c 4 H 2 + c 3 /3 + 12c 4 #) tt 6 
—2 + 3c 4 7r 4 L V / V / 

+36clHf3Tr 7 + 2 (c 2 - 2c 3 tt) - 3c 4 tt 4 (c 2 + 2c 3 7f) - 8Htt (c 3 - 6c 4 tt) + 8c 4 #tt 5 (c 3 + 9c 4 tt)] 

-f 2/3 . . [24c 4 # (l3# 2 + 9H) tt 3 - 108c 2 4 H (2H 2 + 3ij) tt 7 + 2c 2 7f + 3c 2 c 4 7T 4 7f 

— 2 + 3c 4 7T 4 L \ / V / 

+6i^7T (c 2 - 4c 3 7f) - c a Htt 5 (7c 2 + 36c 3 tt) - 12tt 2 (^c 3 i^ + 3# 2 (c 3 - 8c 4 7f)^) 

+6c 4 tt 6 (3c 3 £T + S 2 (5c 3 + 54c 4 7f)) (All) 



= 



